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Abstract 

Local solutions of the multidimensional Navier-Stokes equations for isentropic compress- 
ible flow are constructed with spherically symmetric initial data between a solid core and a 
free boundary connected to a surrounding vacuum state. The viscosity coefficients A, p are 
proportional to p e , < 6 < 7, where p is the density and 7 > 1 is the physical constant of 
polytropic fluid. It is also proved that no vacuum develops between the solid core and the 
free boundary, and the free boundary expands with finite speed. 

1 Introduction 

We study compressible Navier-Stokes equations with degenerate viscosity coefficients in M n (n > 
2), which can be written in Eulerian coordinates as 



d t p + v • (pu) = 0, 

d t (pu) + V-(pu®u)+VP = dW(p(S/u + Vu T )) + V(Adivu) + pf 



where p,P,u = (ui,--- ,u n ) and / = • • • , f n ) are the density, pressure, velocity and the 
given external force, respectively; A = X(p) and p = p(p) are the viscosity coefficients with the 
property A(0) = p(0) = 0. 

In this paper, we focus on the following initial-boundary value problem for (jl.ip . The initial 
conditions are 

(p,u)\t=o = (po,uo)(x) = (po(r),uo(r)-),0 < a < r < 6, (1.2) 



where r = |x| = ^Jx\ + ■ ■ ■ -\- and b > a > are two constants, the boundary conditions are 

U\r=a = 0, (1.3) 

p\ x edn t = 0, (1.4) 

where dQt = ip(dQ,Q,t) is a free boundary. Here, OQq = {x G R n | |x| = b} is the initial free 
boundary and -ip is the flow of u: 

dtip(x, t) = u{tp(x, t),t), x g R n 
^(x,0) =x. 
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For the initial-boundary value problem (jl.l|) - (|1.4p with the spherically symmetric external 
force 

/(x,i) = /(r,i)~ 
r 

we are looking for a spherically symmetric solution (p, u) : 

p(x, t ) = p(r, t ) , u(x, t) =u(r,t)-. 

r 

Then (p, u) is determined by 

d tP + d r {pu) + ^pu = 0, 

Tl — 1 Tl — 1 

p{p t u + ud r u) + d r P = (A + 2p)(d 2 r u H <9 r u — u) q.5) 

77- — 1 

+2d r pd r u + d r X(d r u H u) + p/, 

r 

with the initial data 

(p,u)\ t=0 = (po,uo)(r), < a < r < 6, (1.6) 



the fixed boundary condition 
and the free boundary condition 



U\ r=a = 0, (1.7) 

p(b(t),t)=0, (1.8) 

where 6(0) = b, V(t) = u(b(t),t), t > 0. 

To handle this problem, it is convenient to reduce the problem in Eulerian coordinates (r, t) 
to the problem in Lagrangian coordinates (x, t) moving with the fluid, via the transformation: 

n-1 



x= y n L p(y,t)dy, 

J a 

then the fixed boundary r = a and the free boundary r = b(t) become 

5 = and 6(t) = / y n - 1 p{y 1 t)dy = / y^ 1 p (y)dy, 

J a J a 

where y n ~ l po{y)dy is the total mass initially. Without loss of generality, we can normalized 
it to 1, so that the region {(r, t)\a < r < b(t),0 < t < T} under consideration is transformed 
into the region {(x,t)\0 < x < 1, < t < T}. 

Under the coordinate transformation, the equations (|1.5p ^( TL8|) are transformed into 



d t p = -p 2 d x (r n x n), 

8 t u = r n ~ l d x (p{\ + 2p)d x {r n - 1 u) - P) - 2(n - l)r n ~ 2 ud x p + f(r, t), (1.9) 
r n (x,t) =a n + n p~ l (y,t)dy, 

with the initial data 

(p, u)\ t =o = (po,u )(x),r\ t =o = r (x) = (a n + nj pQ 1 (y)dy^j , (1.10) 
and the boundary conditions: 

u(0,t) = 0, (1.11) 
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P (i,t) = o. 



(1.12) 



It is a well- know fact that if we can solve the problem (|1.9p ^ (|1.12p . then the free boundary 
problem (|l,5p ^ ([1.8p has a solution. 

For simplicity of presentation, we consider only the polytropic gas in this paper. That is, we 
assume 

P{p) = C P \ 1 >l 1 
p(p) = ap 9 , 

x(p) = c 2P 9 ,e>o. 

Without lose of generality, we assume C = 1 . Additionally, we assume that the external force 
/(r, t) satisfies 

f(r,t) £ C l {[a, +00) x [0,+oo)) and \ f(r,t)\ < f(t), V r > a, t > 0, (1.13) 

with / G C x ([0,oo)) and f> 0. 

The results in [U [16] show that the compressible Navier-Stokes system with the constant 
viscosity coefficient has a singularity at vacuum. Considering the modified Navier-Stokes system 
in which the viscosity coefficients depend on the density, Liu, Xin and Yang in [8] proved 
that such system is local well-posed. It is motivated by the physical consideration that in the 
derivation of the Navier-Stokes equations from the Boltzmann equation through the Chapman- 
Enskog expansion to the second order, cf. the viscosity is a function of the temperature. 
If we consider the case of isentropic gas, this dependence is reduced to the dependence on the 
density function. 

In this paper, we establish the existence, uniqueness and the continuous dependence on initial 
data of the local solution to the initial-boundary value problem (]1.9p ^ ()1.12p . The important 
physical consequences of our results are that no vacuum states can occur between the solid 
core and the free boundary, the interface separating the flow and vacuum propagates with finite 
speed. 

The main analytical difficulty to handle this problem is the singularity of the solution near 
the free boundary. The most difficult thing is to get the low bound of p, i.e. p > ^ (1 — x) a . 
Since n > 2 and the viscosity coefficient p depends on p, the nonlinear term r n ~ 2 ud x p in (jl.9p -? 
makes the analysis significantly different from the one-dimensional case [21 [8[ fT5j [T71 [T8] . For 
example, in the proof of Proposition 12.21 Lemmas 12.51 and !5.1[ the nonlinear term do cause much 
trouble. Another difficult thing is to obtain the bound of ||u||^4m. In one-dimensional case, one 
can obtain the bound of ||u||^4m easily, but it's difficult in our case. We apply the inductive 
method to prove it in Lemma [2. 2 1 In the proof of Proposition [2T2l we must get -%rr dy < Cp 9 . 
For this purpose, Yang-Zhao|17| needed to obtain the L norm of dtu in one-dimensional case. 
Here, we only need to get the I? norm of d^u. In this paper, we use the continuation method, 
the energy method and the pointwise estimate techniques to deal with the singularity of the 
solution near the free boundary, and obtain the bounded estimates of the solution. 

Under some new a 'priori estimates, using an effective difference scheme to construct ap- 
proximate solutions, we could obtain the existence of the local solution with < 9 < 7. Using 
weighted function p\p^ (pi — P2) 2 , h + h = — 3, (not only l\ = — 1 + and I2 = —2 as in [3]), 
we prove that the local solution is continuously dependent on the initial data. To do this, we 
need a6 < ^, so 9 can be any value in (0, 7) when a < ^. 

For related free boundary problems for isentropic fluids with density-dependent viscosity (p = 
p 9 ), see Fang-Zhang[2j, Liu-Xin-Yang|8], Vong-Yang-ZhufTp"] . Yang -Zhao[T7] and the references 
cited therein. For the other one-dimensional or spherically symmetric solutions of the Navier- 
Stokes equations for isentropic compressible flow with a free boundary connected to a vacuum 
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state, see MatusuNecasova-Okada-Makino[9j, Nishida jlOl Ill]. Okada[12]. Okada-Makino[13] and 
the references cited therein. Also see Lions |7J and Vaigant[T4"] for multidimensional isentropic 
fluids. Especially, Chen-Kratka[lJ studied the free boundary problem of compressible heat- 
conducting flow with constant viscosity coefficients and spherically symmetry initial data. 

Now we list some assumptions on the initial data and constants (7, 6, c±, C2). 

(Al) Let ci,C2,9, 7 satisfy 

ci > 0, 2ci + nc 2 > 0, 

O<0<7,7>1. (1.14) 
(A2) There exist three positive constants < a < 1, A > B > satisfying 

i-<«<min{A,^l T }, aie-lX 1 -, (1.15) 



and 

(A3) Assume 



A(l-x) a < Po ( x ) <B{l-x) a . (1.16) 
u o eL°°([0,l}), ^+ e (^G^([0,l]), 
(^),€L 2 ([0,1]), (l-xr^GL^O,!]), 



(2c 1+C2 )p^ i (r "-^ )x -2c 1 (n-l)^(^) a; G^([0,l]), (1.17) 

/ x Tq 

where 1 — 2a9 < ao < min{l, 1 + 2a — 2a9}, and m is a integer satisfying 

1 11 

m > max{ , — -, ,2}. (1-18) 

x l + a0-a 2(1 -0a) 4 -4a' s v ; 

Under the above assumptions, our main results can be stated as follows: 

Theorem 1.1 (Existence). Under the assumptions (Al)^ (A3), there is a positive constant 
T\ > such that the free boundary problem $1.9\) ^ ll. ii?J) admits a weak solution (p,u,r)(x,t) 
on [0, 1] x [0, Ti] in the sense that 

p(x,t),u(x,t),r(x,t) €L°°([0,1] x [0,T 1 ])nC 1 ([0,T 1 ];L 2 ([0,l])), 
j (1 - x) a < p(x, t) < 3B(1 - x) a , (x, t) G [0, 1] x [0, Tx], 

p e+1 d x {r n - l u) £L°°([0,1] x [0,T 1 ])n^([0,T 1 ];L 2 ([0,l])), 
^.r,a x «€i 00 ([0,Ti] J i Ao ([0,l])), 
where \q > 1 is a constant which will be defined at i2.19\) . and the following equations hold: 

d t p = -p 2 d x (r n ~ l u), p{x, 0) = p , 
r n (x,t) = a n + n \ p~ 1 (y,t)dy, d t r = u 
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for almost all x G [0, 1], any t G [0, T%\, 




(p)dxdt 



(1.19) 



for any test function 4>(x,t) G C^°([0, Ti) x (0, 1]). Furthermore, if a9 < ^, the solution satisfies 

pd x {r n - l u) G L°°([0, 1] x [0, Ti]), (1.20) 

(l-x)-^GC([0,T 1 ];L oo ([0 J l])). (1.21) 

Theorem 1.2 (Continuous Dependence). Under the assumptions of Theorem \l.l \ and < a9 < 
h, if (pi,Ui,ri) is a solution to the system hl.9\) ^ ll.l2\) with the initial data (poi,uoi,roi), and 
satisfies regularity conditions in Theorem \l.l\ i = 1,2, then we have 

[ ((«i - u 2 f + P \- e pf-\ Pl - P2 f + P W(n ~ dx 

< Ce ct J ((n i-n 02 ) 2 + />orVo2" 4 (/ ) oi-po2) 2 + PoiPo2 1 ^oi-r-o2) 2 )^, 
for all0<t< T 1 . 

Remark 1.1. It is noticed that the set of initial data (po, Uo, tq) verifying all assumptions in 
Theorem 1 1 . 1 1 and 1 1 . 2 1 contains a quite general family of functions. For example, if po = C(l — x) a 
with the exponent a satisfying (jl.15)) . then it satisfies all assumptions on the initial density. 

From Theorem 11.21 we can obtain the uniqueness of the solution immediately. 

Theorem 1.3 (Uniqueness). Under the assumptions of Theorem \1. 6 A if the free boundary prob- 
lem il.9\) ~ [1.12\) has two weak solutions (pi,Ui,ri) (i = 1,2) on [0,1] x [0, T±] as described in 
Theorem \l.ll then for all t G [0, Ti] and almost all x£ (0,1), (pi, u\, r\)(x, t) = (p 2 , u 2 , r 2 )(x, t). 

The rest of this paper is organized as follows. In Section [21 we obtain some basic energy 
estimates, the lower and upper bounds of the density. Then in Section [3j we construct approxi- 
mate solutions by the difference scheme, and obtain the existence of the local solution. And in 
Section [U we prove Theorem 11.21 by the energy method. At last, in Section we prove some 
useful lemmas, which are used in Section [2j In what follows, we always use C(Cj) to denote a 
generic positive constant depending only on the initial data. 

2 Some a priori estimates 

In this section, for simplicity of presentation, we establish certain a priori estimates for smooth 
solutions to the initial boundary value problem (jl.9|) ^ (|1.12p . 
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Proposition 2.1. Let {p,u,r)(x,t), x G [0,1], t G [0,T], be a solution of (fOj) ~ ([77Tg|) . T/ien 
under the assumptions of Theorem \l.l[ there exists a positive constant C\ = Ci(||/Oo||Lg°) H^oIIl 2 ? 
H/IUf ,) such that, for all t G [0,T] 

j {\u 2 + ^—p^ l )[x,t)dx + f f\{ C2 + -c 1 )p e+ \d x {r n - 1 u)f 
Jo 2 7 - 1 Jo Jo n 

+ 2(n ~ ^ cyp 0+1 (r n - x d x u - -) 2 }(x,s)dxds 
n rp 

< e'/ 2 (j\\v?o + ^pV^dx + jT J 2 ( S )d S ) < C ie */ 2 . (2.1) 

Proof. Multiplying (|1.9)) i by p 7 ~ 2 , (II. 9p ? by u and integrating it over [0,1] x [0,t], using the 
boundary conditions (|l.lip ~ (|1.12p and Young's inequality, we have 



+ - 



7- 


1' 


cip e 


-hi 


1 




7 - 


1' 


1 





-t r\ 
, / i 

/o Jo « 



L )(x,t)dx+ / / {( C2 + ^ Cl )p^i(a :c (r- 1 n)) 2 



n r/3 

l )dx + \ \ fudxds 
Jo Jo 

< / (o u + : tpT 1 )^ + 7 / / "" 2 ^s+ / 7 2 c?s. 

Jo 2 7-1 4 Jo Jo Jo 

Using Gronwall's inequality and (Al), we can obtain (|2.ip easily. □ 
Lemma 2.1. Under the assumptions of Proposition we have for all t G [0, T], 

d t r(x,t) = u(x,t), d x r(x,t) = > 0, r > a, (2.2) 

rPdJ{x,t) = (r$d x p e )(x) - -^^((r 1+ ^)(*,i) - (rl + ^ n u )(x)) 

ZCl + C2 

-r^p 7 + r 1+/3 " n / +(! + /?- n)uV" n ) (x, a)ds, (2.3) 



2ci + c 2 Jo 
where f3 = *s a positive constant. 

Proof. From (|1.9|U . we can get (|2.2|) immediately. From (|1.9P i and ([2.2p . we have 

d t {r p d xP e ) = (3rP- l ud xP e +r f, d tx p 9 , (2.4) 
d tP e = -e P e+1 d x {r n - l u). (2.5) 



From (|L9|)r, and (|2T4j) ~ ([23]) . we obtain 

Qt(^d xP °) = Pr^udJ -r^d x {9p e+l d x {r n - l u)) 
= ~—?—{r 1+ P- n {d t u -f) + r?d x p~<). 

2C1 + C2 

Integrating both sides of the above equality in time [0,i], using integration by parts, we can 
obtain (|2.3p immediately. □ 
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In Lemmas I3.3M3.5| we use the continuation method to estimate the lower bound of the 
density. The key point is in the following proposition. 

Proposition 2.2. Under the assumptions of Theorem \l.l\ there exists T\ > 0, such that if 

-p (x) < p(x,t) < 3 Po (x), V t e [0,T'\, (2.6) 

where T' G (0, T{[ is any fixed positive constant, then we have 

^p (x) < p(x,t) < 2p (x), V t E [0,T']. (2.7) 

Remark 2.1. Here T\ can be defined by 

Ti =min{l,Ti,r 2 ,T3,r 4 }, (2.8) 

where T x , T 2 , T 3 and T 4 are defined by (l2T23j) . HjZTM . (I23U|) and respectively. We 

assume that T\ < 1, since we only obtain the local existence of solution. 
From (|2.6|) . we obtain 

a<r <(^a n + n J^{^{l-y) a )- l dy^j n <R<oo, for all t E [0, T'], (2.9) 
where R is a constant. 

Under the assumption (12.6H . we could obtain the following lemmas. The proof of Lemmas 
l2.2N2.6l and Corollary 12.11 can be found in Appendix 5. 

Lemma 2.2. Under the assumptions of Theorem \l.l\ and \2. 6\) . there exists a positive constant 
C2 = 6*2(0, IIpo ||z,£°) II^oIIl**™) II/IIl^) Ci), such that for all t G [0, T'] and k = 1, 2, • • • , 2m, 

f u 2h {x, t)dx + f f p e+1 u 2k - 2 r 2n - 2 (d x ufdxds < C 2 . (2.10) 
Jo Jo Jo 

Lemma 2.3. Under the assumptions of Lemma \2.2\ there exists a positive constant C3 = 
C z {a,R,A,B,\\pl +e {u ) 2 x \\ L ,, WlUr^Cx) such that, for allte[0,T'], 

1 

\u(x,t) -u (x)\ 2 dx <C 3 t. (2.11) 







Corollary 2.1. Under the assumptions of Lemma \2.2\ there exists a positive constant C 4 
C 4 (a, i? ; \\uq\\l^ ,02,0^) such that, for all t £ [0, T'], 



t r l 

1,1,1 .. _ . , 

< C 4 t^=T) (1 _ x )^r. (2.12) 



Jx r 



dxds 



Lemma 2.4. Under the assumptions of Lemma \2.2\ there exists a positive constant C5 
C 5 (a,R, ||(1 -x) a °(p o ) 2 x \\ L i, \\f\\ L ?,C 2 ), such that, forte [0,T'}, 



[\l-xr°(d x p e (x,t)) 2 dx < C 5 , (2.13) 
Jo 



where «o > 1 — 2a9. 
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Lemma 2.5. Under the assumptions of Lemma \2.2\ there exist a positive constant C% which is 
defined by i5.13\) and a positive constant C-j such that, for all t G [0,T"], 



!\d t ufdx + f f 1 \p e+1 r 2n - 2 (d tx u) 2 + p^r^idtu) 2 } dxds < C 6 , (2.14) 
Jo Jo Jo L J 

j (p 1+e {d x u) 2 + P e - 1 u 2 ^dx + j j (d t ufdxds < C 7 (2.15) 
l 

p 6+3 (d x u) 4 dx < C 7 . (2.16) 

Lemma 2.6. Under the assumptions of Lemma \2.2\ there exists a positive constant C% such 
that, for allte [0,T'), 

\\ u \\l°°{[o,i]x[o,t']) < Cs, ( 2 -17) 

I \d x u\ x °dx<C 8 , (2.18) 
Jo 

where Xq is a constant satisfying: 

Remark 2.2. Since m > 4 _ x 4a and a < we have 1 < min{ 4r ^^ fl , ^rjifm } an d the set of Ao 
is not empty. 

Now, we turn to prove Proposition 12.21 
Proof. From (|1.9|) . we have 

2ci + C 2 Jo V r / 2c l + c 2 Jo A V KrU 

' (fr l - n )(y,s)dyds. 



2ci + c 2 y 

Using (pJ2|) . we get 




p e (x,t)+ " [ p\x,s)ds 
2ci + c 2 



o 

> pi(x) dt^-vn -x)ssr - — / —ds 

~ m ' 2ci + c 2 V ; 2 Cl + c 2 J r 

2ci+c 2 J J x P x \r) 2 Cl + c 2 UJULt K ' 

2ci + c 2 7 r 

2 Cl (n - 1)61 f 1 p0[)x / ^ _ ^ _ (2 _ 2Q) 



2ci +c 2 Jo J x 

where m > s^^y and C 9 := 2^C 4 , Ci := 2^a- n+1 | L / | L 
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Using Lemma UM P^]) . ([^6]) and < a < 1, we have 



u(x, s)p (x, s) 
r(x, s) 



ds 



< 



(3pof 



/ \u(x,s)\ds < Cs^^t := C u p e t, (2.21) 
Jo a 



p 9 —^—dxds 
/o Jx r 



< 



f f \d x u\dxds < Cs^^t := C 12 (p ) e t, (2.22) 
Jo Jx a 



and 



n/9 e — ?—dxds 
r 



p 9-^ r -n-l udxds 







< a- n - i (3p ) e ' / \\u\\ L oo\\3(p )- L \\ L ids 

Jo 

< Ct(p ) 9 [ (1 - x)- a dx < C 13 ( Po ) e t. 

Jo 



We define 



)2(2m-l) 



(i - (§) v 



From ([230]) ~ (12T231) . we have 

p°{x,t) + 



2ci + c 2 Jo 



t 2 
p 1 (x,s)ds > (-p ) e , 



(2.23) 



(2.24) 



(2.25) 



for all < t < V < Ti < T x . 

To get a lower bound on p(x,t) from (|2.25p . we need to get a upper bound of the term 
Jq p' y (x, s)ds for sufficiently small t. For this purpose, we set 



Z(t) = / p r (x,s)ds. 



(2.26) 



Similar to the proof of (|2.25p . we can deduce that Z(t) satisfies the following differential inequal- 
ity from (fl~9l : 



Z'{t)) 7 + Z(t) < p[i + (C 9 + Cio + C u + C a2 + C 13 )t^^ (1 - x) e °, (2.27) 



2ci + c 2 
for t £ [0, X"]. Define 



(2 e - l)A 6 



2(2m-l) 



Cg + Cio + Cll + C12 + C13 



(2.28) 



For all i < T' < T 2 , we have 



,o6» 



and 



(Z'(t))- + -^Z(t)<C 14 (l-x 

ZCl + C 2 

p(s,t) <2p (x), ViG [0,1], tG [0,T']. 
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This is the second inequality of (|2.7p . 

Noticing Z(0) = and < 6 < 7, we finally deduce from (|2.27|) that 

> C 14 (l - x) da ((l + -1^-(C 14 (1 - xf^t)^ - 1). (2.29) 



2ci + c 2 ' 2ci + c 2 

By choosing T3 > sufficiently small such that 

Cu (l ~ (l + ^^(^4(1 - x)^)^T 3 ) A ) < ((§)« " (^) A\ (2.30) 
we can get from (l239l ^ (l2T30|) that 

Inserting (|2.3ip into (|2.25|) . we can arrive at 

P° > {\po)\ (2.32) 

provided that < t < T' < Ti. This completes the proof of Proposition 12.21 □ 

Using the pointwise bound for p (|2.6p . we can also derive other estimates for (p, it) in the 
following propositions. 

Proposition 2.3. Under the conditions in Lemma \2.Si we have 

f 

Jo 



\d t p(x,t)\ 2 dx <C, (2.33) 
f o f Q d t (p 6+l d x (r n - l u)) {x, ri) 2 dxdr] <C, t€ [0, T'\. (2.34) 



Proof. From equation (|1.9P i and Holder's inequality, we have 

f \p t (x,t)\ 2 dx = f K^Cr"" 1 !*)) (x,t)f dx 
Jo Jo 

< C I' (p 2 r- 2 u 2 + p A r 2n - 2 {d x u) 2 ) (x, t)d 
Jo 

Using Lemma[231 flUJ), QHft, (|2TT6l) and a < we have 

p 2 r - 2 u 2 (x,t)dx <C u 2 (x,t)dx < C, 
Jo 

rl 

p A r 2n - 2 {d x u) 2 {x,t)dx <C / (p 2 d x u) 2 (x,t)dx 

Jo 

~ C \So p9+3 ( dxU ^ dx + J P^dx) < C- 
From above, we get (|2.33[) immediately. 
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Now, we turn to prove ()2 .34j) : 



t r i 



o Jo 



d t (p e+1 d x {r n - l u)) (x,t) 



dxdr] 



< C 



ff 

lo Jo 



p 20+4 [dx ( r n-l u)] 4 + p 20+2 (^(r"" 1 ^) 2 ) (x, V )dxd V . 



Denote 



t r-i 



D 2 



'0 Jo 
rt rl 



JO 



P 2e+4 (^(r n_1 n)) 4 ) (x, V )dxd7], 
P 29+2 (d tx {r n - l u)f) (x,7 1 )dxdrj. 



From Lemmas E21 [231 §2J)§ and (f2Uj) . we have 

*t rl 



D x < C 

< C 

< c 



p2 e+A ( p -4 r ~4 u 4 + r in-A^Q xU) ^ dxdr} 

(V + p e+3 (d x u)^ dxdr, 
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ft r l 



o Jo 



and 



t rl 



D 2 C 
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ft r l 




{p M+2 (d x (r n l -V)) 2 + p M+A (d x (r n -%u)y)dxdr) 

p 2e + 2 r -4 p -2 u 4 + p 2e + 2 r 2n-A {dxU) 2 u 2\ ^ + ^ 



JO 



< c 

< c. 

This completes the proof. 

Proposition 2.4. Under the assumptions of Lemma \2.Sk we get 

\\[p e+1 d x (r n -\)](x,t)\\ L ~ m]x[0 , TI]) < C, 

dx < C, 



d x (p e+1 d x (r n - l u)) 

I \d xP \dx < C, for all t e [0, T'\. 
Jo 



Proof. From (|1.9j) -? . we have 

d x [p 9+l d x {r n - 1 



u 



d t u 



+ 



dxP 1 

(2ci + c 2 )r™- 1 ' 2ci + c 2 
2ci(n — l)ud x p e 



f 



□ 

(2.35) 
(2.36) 

(2.37) 



(2ci+c 2 )r (2ci + c 2 )r n - 1 ' 
From (|1.13j) and Lemmas I2.5N2.61 using similar arguments as that in (|2.2ip ~ (12.23p . we get 

p^d^r^u) 
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< 



(2ci + c 2 )a n ~ 1 



\dtu\ dx + 



1 



2ci + c 2 



, 2ci(n- 1) / jH , 
^(2c 1 + c 2 ) ^ r "~ 
< C 

Similarly, using (|2.13p and «o < 1> we have 



AM- 



+ 



(2ci + c 2 )a 



n— 1 



|/|dx 



^ |a x (^^(r^u)) I < C + C J \d x {p e )\dx 

< C + C [ (1 - x) a °\(p e ) x \ 2 dx+C [ (l-x)~ a °dx 

Jo Jo 

< C + CC 5 + [ (1 -x)~ ao dx < C. 

Jo 

From (|2.6p . (|2.13|) . ao < 1 + 2a — 2a# and Holder's inequality, we get 



\px\dx < c(J (l-x) ao (p e )ldx 



(1 - x)- ao p 2 - 26 dx 



1 \ 2 

a o +2a-2a0 



< C. 



< C^d-x 

This completes the proof. □ 

Lemma 2.7. Under the assumptions of Lemma \2.2\ we have 

\\r(;t)\\L°°n W iM<C, V t £ [0, T 1 ]. 

Lemma 2.8. Under the assumptions of Lemma \2.2\ and < a6 < \, there exists a constant 
C15 > such that 

\\p9xu\\l°°{[o,i}x[o,t']) < Ci 5 , (2.38) 



\Pq Pt\\L°°([0,l]x[Q,T']) < Cl5 



Proof. From (|1.9j) ? . we get 

„i— n 

2ci + c 2 
2d(n-l) 



^ - (n - l)ur~ n 



p r 



l-n r l 



2ci + c 2 



<9 t u 



2ci + c 2 



r -» u + p -^l-n / p v dx lZ\ dy \ + 



p r 



l-n fl 



2ci + c 2 



-dy. 



(2.39) 



It's easy to show that the first, second and fifth terms of the right side of above equation have 
bounded L°° norm. Similar to the proof of (j2.2ip ~( |2.23p . we get the same result of the fourth 
term. Now we turn to estimate the third term. Using (|2.14p and < a.9 < i, we have 



p Or 1 n f 1 dtu 



2ci + c 2 



„n— 1 



dx 



< Cp-\ \ \d t u\ l dxY(\ - x)2 < C(l - x)2~ ay < C. 
'o 



From above all, we get (|2.38p immediately. From (jl.9P i . (|2.6p and (|2.38p . we can obtain (|2.39p 
immediately. □ 
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3 The proof of Theorem 11.11 

To construct a weak solution to the initial boundary value problem (|1.9p ~ (|1.12p . we apply the 
space-discrete difference scheme method as in [TJ, which can be described as follows. 

Let h be an increment in x, such that Nh = 1 for some N G Z + , and Xj = jh for j G 
{0, • ■ ■ , N}. For each integer N, we construct the following time-dependent functions: 

( Pj {t), Uj {t), rj {t)),j = 0,--- ,N 

that form a discrete approximation to (p, u, r)(xj,t) for j = 0, • • ■ , N. 

First, (pi(t),Uj(t),rj(t)), i = 0, • • • , iV, j = 1, • • • , iV, are determined by the following system 
of 3N + 2 differential equations: 



u i 



Pi = -fiS^m), (3.1) 
= r]~Ha j - 2d(n - V(*>?-i) + /,-, (3.2) 

rf +1 = a" + r i ^A ; (3.3) 
fc=0 Pfc 

with boundary conditions 

u (i) = 0, r (t) = a, 
pi - p]+ e (2 Cl + czWr^ux) + 2d(n - l)^ 11 /^ = 0, (3.4) 

and initial data 

I fjh 1 rjh 

^(°) = r/ Po(y)dy, uj(0) = - / u (y)dy, (3.5) 
" -/(j-l)h " 

where j = 1, • • • , iV, 

p (0) = fn(0), «o(0) = 0, r (0) = a, (3.6) 
r"(0)=a" + n^— -, i = 1, • • • ,JV + 1. (3.7) 

Here, 5 is the operator defined by Svjj = (wj+i — Wj)/h, and 

aj(t) = (2 Cl + c 2 )^i 1 1 5(r;r 1 1 n J _ 1 ) - p]_ v fj(t) = f(r j: t). (3.8) 

The basic theory of ordinary differential equations guarantees the local existence and unique- 
ness of smooth solutions (pi,Ui,Vi), i = 0, • • • , N, to the system ()3.1j) - ()3.4[) on some interval 
[0,T h ], such that 

< pi(t) < 00, \uj(t)\, \rj(t)\ < 00, t G [0,T h ], 
where i = 0, • • ■ , iV and j = 0, • • • ,N + 1. Here, we use the fact min pdO) > 0. 

i=0,...,N 

So, we can let T^ ax be the maximal time such that the smooth solutions exist on [0, T^ a:E ) 
and satisfy 

< Pi {t) < oo, |«j(t)|, Ir,-(t)| < 00, t G [0,2^), 

where i = 0, • • ■ , N and j = 0, • • • , N + 1. Our first goal is to show that T^ ax > T\, and the 
solutions satisfy 

»Pi{0) < Pi (t) <3 Pi {0), i = 0,---,N, (3.9) 
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for all t € [0, Ti], where T\ > is given in (|2.8p and independent of h. 

Based on the work of Proposition 12.11 and Lemma 12.11 using similar arguments as that in 
[HE], we have following two lemmas. 

Lemma 3.1. Let (pj{t),Uj(t),rj(t)), j = 0, • • • ,N be the solutions to (JX7J) ~ Qff.^p , then there 
exists a positive constant C(A, ||ito||z,2 , /) such that, for all t 6 [0,T^ ax ), 

j=o v ' J0 j=0 

2(n - 1 



/ 1 i \ / 

e (?'? + * + / E (<** + «) (* (-r 1 - 

pJ+'OT'fc'j - ^«) 2 ) * < Ce'. (3.10) 



Lemma 3.2. XYie solutions (pi(t), Ui(t), rj(t)), i = 0, ■ • • , JV, satisfy the following identities: 



n 



d tri (t)= Ul (t),6r?(t) = - 

P: 



rpp%t) = r? (0)^(0) - (rj+^^C*) - rj^-(0)^(0) 

f jT (-rf 5pJ + + (! + /?- n)u?rf- n )( S )<fe 



2ci + c 2 

for all i = 0, ■ ■ ■ , N , j = 0, ■ ■ ■ ,N — 1 and t E [0, T^ a J. 

Based on the work of Proposition 12.21 and 12.41 using similar arguments as that in [HE], we 
have following lemma. 

Lemma 3.3. Under the assumptions in Theorem li.il we have that for all h € (0, ho] (ho is a 

sufficiently small positive constant ), there is a T\ > such that, if 

\pm<Pi{t)<z P m, 

for all i = 0, . . . , N and t G [0, T'] where T' € (0, T^ ax ) n (0, Ti], then we have 

< ipi(0) < Pi (t) < 2(h{0), (3.11) 
a<rj(x,t) <C, (3.12) 
\uj(x,t)\<C, (3.13) 

for all i = 0, • • • , N, j = 0, ■ • • , N + 1 and t € [0, T'], where C is independent of T' and h. 

Remark 3.1. For simplicity of presentation, in Section [21 we establish some a priori estimates 
in the continuous version to the initial boundary value problem (jl.9p ~ (|1.12p . so we need ho is 
a sufficiently small positive constant. 

Then, applying the continuation method, we obtain the following lemma. 
Lemma 3.4. Under the assumptions in Lemma Wlh we have i3. 9\) holds for all t G [0,T^ ax ) n 

[o.vV- 



14 



Proof. Let A = {T 1 G [0, T^ ax ) n [0, 2\] | (02} holds for all t G [0, T'] }. 

Since minyOj(O) > 0, we have G C([0, T^ ax )), i = 0,...,N, and there exists G 

(0, T^ ax ) such that (|3.9p holds for all t G [0, T^]. Thus, A is not empty and relatively closed in 
[0, T^ ax ) n [0, T\\. To show that A is also relatively open in [0,T^ ax ) n [0,Ti], and hence the 
entire interval, it therefore suffices to show that the weaker bound 

^Pi(0) < Pi (t) < 3 Pi (0), i = 0,...,iV, i G [0,T'] C [0,^)1-1 [0,Ti], 

implies ([3TTTD holds for all t G [0, T'] . Thus, from LemmaEl we have A = [0, T^ ax ) n [0, Ti] . □ 

Then, applying the reduction to absurdity, we have the following lemma. 
Lemma 3.5. Under the assumptions in Lemma VS.'A we have T^ ax > T\. 

Proof. If T^ ax < Ti, from Lemmas E3H331 we have (^TT]) ~ (|3TT3"|1 hold on t G [0,2^). 
Thus, we can extend the existence interval [0,T^ ax ) to [0,T^ ax ]. It contradicts the definition of 
T h □ 

- 1 max • 1 — 1 

Based on the work of Lemmas l2.2N2.8l using similar arguments as that in [HE], we have 
following lemma. 

Lemma 3.6. Under the assumptions in Lemma \3.3[ we have that (pj(t),Uj(t),rj(t)), j = 
0, • • • ,N, satisfy 

N t N 

1 £uf n (t)h + / £ [p e +\s)uf l - 2 {s)8{r]-\s)u j {s)f 

3=0 J ° j=0 

+p e - l (s)rj 2 {s)u* m {s)\ hds < C, (3.14) 

N-l N N 

E (1 - jhr{5p%t)fh + "£(d tUj ) 2 h + E l^-^l^/i < C, (3.15) 
i=o j=o i=o 

/icft < C, (3.16) 



i=o J u i=0 

llft ?+1 ''('-rV)(t)||^([o,T 1 ]) + K(i)|U°o ([0)Tl]) < C, (3.17) 



iV-1 N-l 

E i^r^r 1 ^)) w + E ^ c ' ( 3 - 18 ) 

j=0 j=0 



d 



IriW+YslSr^h+l-nlKC, (3.19) 
i=o 

for i = 0, • ■ ■ , TV, < s < f < T\. Furthermore, if a6 < ^, we /iawe 

ll«%r 1 ^)llL-([o,r 1 ]) + lk- 1 (0)^|U«([o^]) i = 0,---,N. (3.20) 
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Now, we turn to prove Theorem 11.11 

With (pi, Uj, rj), i = 0, • • ■ , N, we can define our approximate solutions (p N ,u N , r )(x, t) for 
the system (jl.9p ~ (jl.l2j) . For each fixed N and i € [0, T\], we define piecewise linear continuous 
functions (p N ,u N ,r N )(x,t) with respect to x as follows: when x G ([xN]h, ([xN] + 



= P[xN] (t) + (sJV - [xN])(p [xN]+1 (t) - P[xN] {t)), (3.21) 
u^a;,*) = u [xN] (t) + (xiV - [xN])(u [xN]+1 (t) - u [xN] (t)), (3.22) 

»t / \ 1/n 

r"(3,t) = (rf xN] (t) + (xN - [xN])(rf xN]+1 (t) - rf xN] (t))) . (3.23) 
We have for jh < x < (j + l)h 

d x u N (x,t) = Uj+l{t) - Uj(t) . (3.24) 
h 

We also introduce the corresponding step functions: 

(Ph,Uh,r h )(x,t) = (p[ xN ],U[ xN ],r [xN] )(t), x G ([xN]h, ([xN] + l)h}. (3.25) 



Using Lemmas I3.lti3.61 we have, for all h E (0, /io] 

SUp {\\u N (;t)\\ LoanWl ,, + ||«f (.,t)|| £2 + \\d t ( Ph )(;t)\\ L2 ) < C, 

te[o,n] 

^ P N (x,o) < P N (x,t) < 3 P N (x,o), v (x,t) e [0,1] x [o,r x ], 
sup f / 1 (i-xr(( / , iV ) e )^ + rv' [0)1] (^) N ) <c, 

te[o,Ti] wo / 

sup (||(p^a«((r^) n - 1 u^))|| L? +rV[o 1 i](pi +fl a«((r JV ) n - 1 ^))) <c, 
te[o,Ti] v 7 

/ Tl ||ft(pl +fl « B ((r JV ) w - 1 u JV ))(. I t)||i a ctt < C, 
io 

sup (^^(^^K^n^i.Ao + ||9 t r iV (-,t)|| L cx 1 ) < C, 
te[o,Ti] 

where TV(g) is the total variation of g, x G [0, 1] and t G [0, T\\. 

With the above estimates, using Helly's theorem and similar arguments as that in Section 9 
of pQ , we can get the following compactness of the approximate solutions (If necessary, we can 
choose the subsequence.): 

(Ph,p N ,u N , r N )(x, t) -> (p, p, u, r)(x, t), a.e. 

{p h ) l+e d x {(r N ) n - l u N ){x,t) - p 1+e d x (r n - x u){x,t), a.e. 

(d x r N ,d x u N )(x,t) -± (d x r,d x u){x,t), weakly in L A °([0,1] x [0,Ti]), 

d t p N (x, t) d t p(x, t), weakly in L 2 ([0, 1] x [0, T{\), (3.26) 

when iV — > co, where r(x, t) = (a n + n J*^ ^^y)"- 
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For any given test function ip(x,t) € Cg°((0, 1) x [0, Ti]), we choose N = ^ is large enough 
such that the support of the test function tp is away enough from the boundaries, that is, supp 
V'C((/ l ,l-/ i )x[0,T 1 ]) = ((i,l-i)x[0,T 1 ]). ' 

Define 

^(f) = ij} h {x,t) = ip([xN]h,t), ih<x<(i + l)h. (3.27) 

We can see that ^(t) = for i = 0,N - 1, N. 

Multiplying (|3.1|) by ipi, summing it up for i = 0, • • • , N, and integrating it over [0, Ti], we 

get 

Ti N 

{i>idtPi + ^iplSir^Ui)) hdt 

i=0 

= 0(h) + ^ ty h d t p h + MPh) 2 d x ((r N ) n - l u N )) dxdt 
Jo Jo 

= 0(h)+ r [ {^d tPh + ^{p h fd x {{r N ) n - l u N ))dxdt. (3.28) 
Jo Jo 

Using (|3.26p . we can pass the limit in (|3.28j) to obtain the first equation in (jl.9j) . 

For any given test function 4>(x,t) G Cg°((0, 1] x [0,Ti)), we choose N = is large enough 
such that the support of the test function <p is away enough from the fixed boundary, that is, 

4>\xe[o,h] = 0. 
Define 

<j>i(t) = 4> h (x,t) = <j>([xN]h,t), (3.29) 

<f> N {x,t) = <P [xN] (t) + (xN - [xN])(<t> [xN]+1 {t) - cj, {xN] (t)), (3.30) 

when ih < x < (i + l)h. We can see that <j)i(i) = for i = 0, N — 1, N. 

Multiplying (|3.2p by (j)j, summing it up for j = 1, ■ ■ ■ ,N, and integrating it over [0, Ti], we 

get 

, Tl N „ Tl N 



= / 1 J2(j) j dtUjhdt- f * y^fyr^Scjjhdt 
Jo j=1 Jo j=1 

Tl N Tl N 

+2ci(n-l) / L,<i>jr?- 2 u j 6(p j _ 1 )hdb- / J2&fi hdt 
Jo =1 Jo - =1 



3=1 " U 3=1 

0(h)- [ cf) h u N (x,0)dx - [ 1 f d t <p h u N (x,t)dxdt 
Jo Jo Jo 



+ 



[ Tl [\2c 1 + c 2 )d x ( ( f> N r N )p e h +1 d x ((r N r- 1 u N )dxdt 
Jo Jo 

Tl f d x ( ( j ) N r N )pldxdt- f^ctin -^(d^i^T- 2 ^ 
o Jo Jo Jo 

+ (n - 2)<p h (r N ) n - 3 d x r N u N + ^{r^ 2 d x u N ) p{dxdt 

/ 4>hfhdxdt, (3.31) 
o Jo 

where fh(x,t) = f\ x m(t), x G ([xN]h, ([xN] + l)h]. Using (|3.26p . we can pass the limit in (|3.3ip 
to obtain (|1.19p . Thus we complete the proof of Theorem 11.11 



17 



Remark 3.2. When aO <\, from ([ODD , we have 

||p Ar ((r Ar )"- 1 n JV ) :[ .|| L o O([0il]><[0iTl]) + \\{ P N (xM~\p N )t\\L^i Wl] ) < C. 
Thus, the limit function (p, u, r) satisfies 

||/j(r n_1 u) x || L oo([ 0il ] x [ 0)Tl ]) + ||pg Pt |U°o ([0,1] x [0,71]) ^ Cj 

and 

(l-xrVGCQO,^];^^,!])). 
Remark 3.3. There is another method to prove Theorem ll.il At first, we consider system 

dtp = —p 2 d x (r n ~ 1 u) 

d t u = r n ~ l d x (p{\ + 2p)d x {r n - 1 u) - P) - 2(n - l)r n ~ 2 ud x p + f(r(x, t),t) (3.32) 
r n (x,t) = a n + n J* p~ 1 (y,t)dy 

with the initial data 

(p,u)\ t=0 = (p £ ,u £ )(x),r\ t=0 = r[ ] (x) = ^a n + n^ ( P l)~\y)dyY , (3.33) 

where p§(> £)>'"o' r 6 converge to po,uo,ro in some suitable spaces as e goes to 0, and the 
boundary conditions: 

u(0,t) = 0, (3.34) 

{p(A + 2p)d x {r n ~ l u) - P - 2(n - V} | x=1 = 0. (3.35) 

Using similar arguments as that in Sections [2H21 we can obtain the existence of the weak 
solution (p £ ,u e ,r e ) to the system (|3.32[) ~ (l3.35f) . and some uniform estimates of the solution. 
Letting e — > 0, we can prove that the limit function (p, u, r) is the weak solution to the system 

ra~ (fTT2i . 

4 Continuous dependence 

In this section, we will prove Theorem ll.2l applying the energy method. Let (pi, ui, r\){x, t) and 
(p2, U2, r 2 )(x, t) be two solutions in Theorem [TTT] corresponding to the initial data (poi, uqi, r i)(x, i) 
and (po2, ^02, ^02)(a ; 5 i) 5 respectively. Then we have, i = 1,2, (x,t) € [0,1] x [0, T{\, 

A 

-(l-x) a <Pi(x,t) <3B(l-x) a , \ui{x,t)\<C, a<n<C. (4.1) 
From Lemma l2.8t we can easily get 

llp^rf-V)!^ ([o,i] x [0,21]) < C. (4.2) 

For simplicity, we may assume that (pi, u%, ri)(x, t) and (p2, f^X^j i) are suitably smooth 
since the following estimates are valid for the solutions with the regularities indicated in Theorem 
11.11 by using the Friedrichs mollifier. 

Let 

q = pi - p 2 , w = u 1 -u 2 , 71 = ri - r 2 . 
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£0 = P01 ~ P02, W = Uqi - UQ2, ^-0 = ^01 ~ ?*02- 

From (|1.9|) i . (|4.ip ^ (|4.2p . Lemma 12. II and Young's inequality, we have 

= [ 2p e 1 p 2 1 KK t - 9p e l +1 d x (r^ 1 u 1 )p 2 1 K 2 + p\d x {rl- l u 2 )Tl 2 dx 
Jo 

< C £ [ p\p 2 x 1l 2 dx + E [ ^jp^dx, 
Jo Jo r l 



where e > is chosen later. 
From (|1.9ft i . we have 



d te = -plr?- 1 d x w- 6 (d x u 2 rr 1 (pi + P2) + (n-l)—) 

-(n -l)^w- pld x u 2 {rr l ~ rT 1 ) + {n- 1)^±K, 
r\ ~ r±r 2 



and 



d 
df 



f p\- e pf^ e 2 dx 
Jo 

f r, 1-0 26»-4 o / 1 Q \ 2-Qp, / n-1 n 26»-4 2 

= / 2 Pi P 2 QOtQ ~ (1 - 0)Px o x {r 1 ui)p 2 q 
Jo 

-{29 - A)p\- e d x {r^- l u 2 )p 2 2 e ^ 6 2 dx 

< 2 f 1 p^pf^edtQdx + C f 1 p\- e p 2 2 e - 4 Q 2 dx 

Jo Jo 

= 2 f 1 -pl-°pf-*grr l d x w - p\- e pf^Q 2 (d x u 2 rr l (Pi + fb) + (n - 1)^) 

Jo ' r 2 

-in - i f 1 ' P r f 6 w - p\~ d p 2 t 2 e d x u 2 {rr l ~ rT 1 ) 

+{n-l) U2pl 92 e Kdx + C p\- e p 2e - A 6 2 dx 
nr 2 Jo 

< C £ f 1 p\- e p 2 2 e - A Q 2 dx + C £ f 1 p^p e 2 K 2 dx + e f'^p^dx 

Jo Jo Jo r l 

+e I'' p^rf^^wfdx. 
Jo 

Remark 4.1. We can use the weighted function p£ p l 2 (pi — p 2 ) 2 , where l±, l 2 satisfy l\ + l 2 
For simple, we choose l\ = 1 — 9, l 2 = 29 — 4. 

We only give a part of the proof of the last inequality. The rests are the same. 

f -p\- Q p 2 2 e ~ A Qrl- l d x wdx < BJ 1 - 1 f p^d x w . pf ~p e 2 - 2 Q . p\- e p e 2 - 2 dx 
Jo Jo 

or e+i i-e 

< R n -\—) 2 - e J W d xW . Pl > p e 2 - 2 e \dx 
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< a f 1 p\ +1 r\ n - 2 (d x w) 2 dx + C e r p\- e pf- 4 g 2 dx. 
Jo Jo 

From equation (|1.9j) '? and boundary conditions (|1.10p - (|1.12p . we get 

1 )u?(x,t)dx = ! {-(\(p 1 ) + 2p(p 1 ))p 1 d x (r^ 1 u 1 )d x (r^ 1 w) 



dt L 2 



+(A(p 2 ) + 2p(p 2 ))p 2 d x (r^ 1 u 2 )d x (r^- 1 w)} dx 



+ / {P(p 1 )d x {r™- l w)-P{p 2 )d(r™- 1 w)}dx 



+2(n-l) / {p( Pl )d x (r^- 2 u lW )-p(p 2 )d x (r^- 2 u 2 w)}dx 
Jo 

+ I (/(n,t) - f(r 2 ,t))wdx := U! + U 2 + U 3 + U A . (4.6) 
J o 



Using the similar argument as that in Proposition ^. 11 Cauchy-Schwartz inequality and (|4.ip ~ (|4.2p . 
we have 

U x + U 3 < -~ f {{c 2 + ~c x )p\ +6 [d x {r«- l w)] 2 
1 Jo n 

+ 2 Cl (n-l) [r ni dxW _^ ?}dx 

n T\p\ 

+c !\p\- 6 P f- A e 2 + P WK 2 )d x , 

Jo 

U 2 < e f p\ +e r 2n - 2 [d x w] 2 dx + e f p\- l w 2 dx 
Jo Jo 

Jo 



(4.7) 



(4.8) 



and 



U 4 = f (f(n,t)- f(r 2 ,t))wdx (4.9) 
Jo 

< C \Kw\dx<C £ l p[' l n 2 dx + e I p\~ e w 2 dx 

Jo Jo Jo 

< C £ [ pf-^dx + eCWwWlco 

Jo 

< C £ f 1 pl-Wdx + eC f 1 p{ +1 r 2n - 2 {d x wfdx f p- x e - x rf n + 2 dx 

Jo Jo Jo 

< C £ f p\- l 1l 2 dx + eC f p\ +1 rl n ~ 2 (d x w) 2 dx. (4.10) 

Jo Jo 

Here, we use the fact that f(r, t) € C 1 ([a, +oo) x [0, +oo)). Choosing a sufficiently small positive 
constant e, we obtain 



j^( W 2 +p\- p 2 2 e - 4 e 2 +pWK 2 ) dx+C l6 J\\ ( Pl r^- 2 (d^f + ^j dx 



< C 17 £ (w 2 + pl- e pf-* 6 2 + pWn 2 ) dx, 
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where C%§, Cyj are two positive constants dependent on e. 
Using Gronwall's inequality, we have for any t G [0,Ti], 

[ (w* +P l-°p?-* g * +P WK 2 ) dx < Ce™^ (wl+p^pf^Ql+plp^nt) dx. 
Then, we finish the proof of Theorem 11,21 

5 Appendix 

5.1 Proof of Lemma 12.21 

Proof. We apply the inductive method to prove this lemma. 

First we consider the case of k = 1. From Proposition 12.11 we obtain (|2.10p with k = 1 
immediately. 

Assume (|2.10j) holds for k = I — 1, i.e., 

f 1 u ^) d x+ f f 1 p e+1 u 2 ^- 2 r 2n - 2 {d x u) 2 dxds<C. 
Jo Jo Jo 

Thus, using a < j^g, we have 

||u ,-1 |||oocfa < J (J Idviv 1 - 1 )]^ ds 



< 



j n p e+i (d^u 1 - 1 )) 2 cb\ n P - 6 - i dx^j ds < c. (5.i) 



Now we need to prove (|2.10p holds for k = I. Multiplying (|1.9|) '? by u 21 1 and integrating it 
over x from to 1, using the boundary conditions (jl.lip ~ (|1.12|) . we have 

^u 2l (x,t)dx = -(2c 1 + c 2 )^ 1 / +1 ^(r"- 1 n)9 :c (r"- 1 ^- 1 ) ( ix 



+ f p 1 d x {r n - l u 2l - 1 )dx + 2c 1 {n-l) f p 6 >d x (r n - 2 u 2l )da 
Jo Jo 



+ / fu*- l dx 





Set 
Thus 



:= Gi + G 2 + G 3 + G 4 . (5.2) 
A 2 = p e+l u 2l - 2 r 2n - 2 (d x u) 2 > 0, A 2 = p e - l r- 2 u 21 > 0. 



-(21 - l)(2ci + c 2 ) I A\dx - 2(n - l)Zc 2 f A x A 2 dx 
Jo Jo 

+ (2(n-l)(n-2)ci-(n-l) 2 (2ci + c 2 )) / A\dx. 

Jo 



Then, from (|5.2p and Young's inequality, we get 



d 
dt 



[ ^-u 2l dx + (21 - l)(2ci + c 2 ) / A\dx 
Jo 2t Jo 
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< ef A\dx + C e [ A 2 2 dx+ [ p 1 d x {r n - l u 2l - l )dx+ [ fu 2l ~ l dx. 
Jo Jo Jo Jo 

From (|1.15p . (|2.6p . (|5.ip . Young's inequality and r > a, we obtain 

f A\dx < C C(l-x) a ^a- 2 u 2l dx 
Jo Jo 

< C [\l-x) a ^u 2l - 2 dx + C f'u^dx 

Jo Jo 

< C||u 2 '- 2 || L oo f (l-x) Q(e - 1) dx + C||u 2 || L oo f u 2l dx 

Jo Jo 

< Cdlu^Hloo + llullloc / u 2l dx), 



(5.3) 







f p^d x (r n - l u 2l - l )dx < e I A\dx + C £ + C £ [ 
Jo Jo Jo 



u 2l dx, 



and 



i r i 

fu 2l - l dx <C + C j u 2l dx. 

Jo 



Inserting above inequalities to (|5.3p and choosing e = ( 2; i)(2ci+c 2 ) i ^ wg that there exist 
two positive constants C\% and Cig such that 



4- r U 2l dx+ I' 1 p e+1 u 2l - 2 r 2n - 2 (d x u) 2 dx 
dt Jo Jo 

< Cis + Ci 9 (||^ 1 |||c + ||«|||oc / u 2l dx). 







Then, using Gronwall's inequality, (A3) and (j5,l|) . we get 

u 2k (x, t)dx + J* J (p e+1 u 2k - 2 r 2n - 2 (d x u) 2 ^J (x, s)dxds 

< (C 18 + Ci 9 / \\u l ' 1 \\ 2 L oods)e Cl9 ^MU ds <C 2 , < t < T'. 
Jo 

So (|2.10p holds for k = I. From all above, we get the result immediately. □ 

5.2 Proof of Lemma f2T3l 

Proof. From (jl.9p . we get 

d t (u - «o) = r n - x d x ((2 Cl + c^+^r^u) - p?) - 2c x {n - \)r n - 2 udj + /. 

Multiply the above equation by u — uq and integrating it over x from to 1, similar to the proof 
of proposition l2.lt and using the inequality ab < + ^b q , where - + - = l,p, q > 1, a, b > 0, 
we have 

Jtl V u - U ^ dx 

+ [\c2 + - Cl )p e+1 (d x (r n -\)) 2 + ^-^- Cl p 6+1 (r^d x u --) 2 dx 
Jo n ' n V r PJ 
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(2 Cl + c 2 y +1 d x (r n - l u )d x (r n - l u)dx 

2ci(n - 1) ((n - 2)r~ 2 p~ 1 uu + r n ~ 2 d x uu + r n ~ 2 u(u ) x ) p 9 dx 

"1 f'\ f'\ 

d^-^p^dx- / d x {r n ~ l u Q )p 1 dx+ / (u-u )fdx 
) Jo Jo 

< 1 f 1 p^(d x (r n - 1 u)) 2 dx + - f 1 p ^u 2 r- 2 dx+- [' p e+1 (d x u) 2 r 2n ' 2 dx 

+C e [ p e+l {d x {r n - l u )f dx + cj p e - l u 2 r- 2 dx 
Jo Jo 

+C e I' p e ~ l r- 2 u 2 dx + C e I' p e+1 r 2n - 2 ((u ) x ) 2 dx + C t C p 2 ^ e+1 Ux 
Jo Jo Jo 

+\ I P 9+1 [Ur^uo)) 2 dx + \ I p 2 ^- (e+1 Ux+ f \(u-u )f\dx. 
1 Jo 1 Jo Jo 

Choosing e = min{c 2 + f ci, ^^^ci}, using n > 2, (fLT5|) . (USD and p£ +fl (M )x G Ll « we * 

4 ( l \{u-u fdx+ [\ C2 + -c 1 -e)p e+1 {d x (r n - l u)) 2 dx 
at J Q l Jq n 



+ (\^^ Cl - e)p e+1 (r n - l d x u - - ) 

Jo n V r p; 



< c e f 1 p^id^-^fdx + a f 1 p< 

Jo Jo 



e - l ulr- 2 dx 



-C t f p e+l r 2n - 2 (u )ldx + C e f p 2 ^ 6+1 Ux+ f \(u-u )f\dx 
Jo Jo Jo 



< C. 

Thus, we can get (|2.1ip immediately. 
Now we turn to prove Corollary I2.lt 

Proof. Using integration by parts, we get 

rt rl d t u 




n-l 



dxds 



< 



Jx r 

1 \U - Uq\ 



„n— 1 



dx 



r n—l 



dx + 



1 Uq rt rl 

— j dx + (n — l)u 2 r~ n dxds 

r Jo Jx 

dx + (n — l)u 2 r~ n dxds. 

Jo Jx 



Then, using Lemmas l2.2N2.3l and Holder's inequality, we have 



{n-l) [ [ u 2 r~ n dxds < (n - l)a~ n f [ u 2 dxds 

Jo Jx Jo Jx 

rt rl 

,4mj„\5i-. 



2m-l 



< {n-l)a- n /(/ u Am dx)^dt{\ -x)tst < Ct(l -x)^r, 
Jo Jo 



1 Hv ^- U f l dy < f H y,t) - , (y)| 2 ^)-(l - xf^ 

r n L (y,t) Jo 
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< a l - n { \u{y,t)-u Q {y)\ 2 dy)^=T) 

Jo 

f a m—1 2m-] 

x( / \u(y,t) - u (y)\ dy) 2m( - 2 ™-v (1 - x)^~ 
Jo 

1 , , 2m-l 

< Ct 2 ^™- 1 ) (1 - X)~5ZT, 



and 



' r o 



dx 



< ||no||L-a 2 - 2 " f [r^-r^da 

J X 



< 



L°°CL 



2-2n 



dtr n ~ l ds 



dx 



1 rt 



2-2n D n-2 



< (n-l)||uo||L-a^ n i?' 



i Jo 



t rl 



< C\\u Q \\ L ™a z - Zn R n - z (/ n 4m dx)4^(l-x)^ds<Ct(l-x) 

Jo 

Finally, we can get (|2.12p immediately. 
5.3 Proof of Lemma 12.41 



□ 



Proof. Multiplying (|2.3p by (1 — x) ao r^d x p 9 and integrating it over x from to 1, using Young's 
inequality, we have 
•l 

(1 -x) Q0 (r^p e ) 2 dx 



< 



I f\l-x) ao (r^d x p e ) 2 dx + C f (l-x) ao (r$(p e ) x fdx 



+C {l-x) ao {r l+l3 ' n ufdx + C {l-x) a °{rl +p - n uvfdx 



+c / (i-x) c 



(j\-rPd xP i + r 1+/3 - n f +(!+/?- n)u 2 r p - n )d^j dx. 



From Lemma [221 ao > 1 — 29a and m > 2 (i-0a) ' we nave 



and 



/ (1 - x) QO (r 1+/3 -%) 2 dx < C / (1 - x)^r ao dx + C T u 4m <ix < C, 
Jo Jo Jo 



(l-x) a °( uV _n (l + /3-n)ds) 



< C / (l-x) a °( u 2 dsYdx<C (l-x) ao / u 4 dsds 



l 



< C (l-x)— ao dx + C / u Am dsdx<C. 



JO 



Thus, we get 



(l-x) a °(r l3 d x p e ) 2 dx<C + C (l-x) a ° \r^d x p i \ 2 {x,s)dsdx 
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< C + C f max(p 7 -") 2 f (1 - x) a ° (r l3 d x p e ) 2 dxds 
Jo Jo 



t r l 



< C + C 



{l-x^^d^fdxds. 



Jo 



Using Gronwall's inequality, we can obtain (|2. 13|) immediately.. 



□ 



5.4 Proof of lemma T2.5I 

At first, we consider the following two lemmas. 

Lemma 5.1. Under the assumptions of Lemma \2.SX we have, for any t 6 [0, T'\, 



t r i 



[p e - l u 2 + p e+l {d x u) 2 ^ dx + I J (d t u) 2 (x,s)dxds 



t r l 



< c 



p 3+e (d x u) A dxds. 



(5.4) 



Jo 



Proof. Multiplying (|1 .9[) -? by dfU, integrating it over [0,1] x [0,t], and using the boundary con- 
ditions (fTTT|) - (fl~T2ll . we have 



n(dtu) 2 (x, s)dxds 



t rl 




p 1 d x {r n (d t u))dxds 



o Jo 



t rl 




'0 JO 

rt r l 



+ 

d 




(2ci + c 2 )p L+y d x (r n - L u)d x (r n ~ L (d t u))dxds 

rt r l 

2ci(n - l)p e d x {r n - 2 u{d t u))dxds 



JO 




f(dtu)dxds 



o Jo 



(5.5) 



i=i 



Using (A3) , (|2.ip , (|2.6p , (|2.9p , Cauchy-Schwarz inequality, Young inequality and m > 4 ^ 1+ ^g_ a ^ , 
we have 



H 2 + H 3 
i 



o 



2 -^P^p 1+9 [d x {r n - l u)} 2 + d(n - lJ^Cr^V) 



dx 



+ 



r n - 2 u 2 ) 



[ [ {(2c l +c 2 ){n-l)p l+e d x {r n - l u)d x { 
Jo Jo 

- (2C1 +C2) (1 + ^V+U^u)] 3 + 2^ Cl (n - lj/^fi^V* 
2 r 

-6cw(n - l)p e ^d x {r n - l u) + 2nci(n - l)(n - 2)p e ~ 1 ^ 



u 



-3ci(n - l)(n - 2) / o y ^a r (r"- i u)} f fe ( fe 



< C + C 



J 
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dxds 



-c- 1 



< c + c 



p l+e {d x uf + p e ~W 



dx 



t r l 



JO 
1 



4ra(«-l) 

1i 4 ™ +p T^r +p 3 + e (d x u) 



dxds 



o 



< C + C f I p 3+e (d x u) 4 dxds 
Jo Jo 
-1 



-C 



-i 



p l+0(^)2 + p 0-l u 2 



dx, 



(5.6) 



Hi 



p 1 d x {r n L u)dx 



t fi 





+ / / 7P 7+1 [^(r n_1 u)] 2 cixds 

is Jo Jo 

t rl 



JO 



2{n-l)pi-d x {r n - l u)dxds+ / / n(n - I) p^^ dxds 



JO 



< / p 7 ^(r n - i u)ci2; + C, 



1 



H±<- I / (d t uydxds + C. 
2 Jo Jo 



Using (|5.5p - (|5.8p . we can obtain (|5,4p immediately. 



(5.7) 

(5.8) 
□ 



Lemma 5.2. Under the assumptions of Lemma \2.2\ there exist two constants C20 and Ci\, such 
that, for any t S [0, T'] 



P* +6 (d x u) 



l dx < C 20 + C21 (J (p l+6 {d x uf + V + (d t uf) dx 
J\pd x ufdx < C20 + C21 (J [p 1+0 {d x uf + /~V + {d t uf) dx 
Proof. From (|1 .9^ *? and (|1.12p . we have 

i- 1 d t u 



(5.9) 
(5.10) 



d r u 



p<- 9 - 1 - (n - l)p^ ur' n 



p r 
2ci + C2 



2ci + C2 

2ci +c 2 \pr n J x r 



+ 



2ci + c 2 7a; r n_ 

Using r > a, 7 > 1 and 7 > 0, we get 
-1 

/ +3 (d^) 4 cfo 



T dy. 



(5.11) 



< 



C + C [ p e - 1 dx + C [ p e - l u A dx + C [ p- 39 - 1 ( [ \d t u\dy) dx 

Jo Jo Jo \Jx J 
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+C^p- 3e - 1 ^y ^dyul + lp^uldy] dx+C J P' 39 ' 1 (J dx - 
Then, using a(9 — 1) > —1, a < min{g^, 3^W} and m > 1+Q g_ a ; we have 

f 1 p 6 - l dx < C [\l - x) a{e - l) dx < C, 
Jo Jo 

f /"Vtfa <C I p {e ' l) ^dx + C I u Am dx < C, 
Jo Jo Jo 

£ p- 36 - 1 \d t u\d^j dx<C j\-™- x (^j\dtu\ 2 dy)*{l-x)£\ dx 
< c(j\d t u\ 2 dx\ j\l-x)~ a ^ 6+1 ^ 2 dx<c(j\d t u\ 2 dx\ , 

jf 1 p- 39 - 1 (T \f\dy\ 4 dx < C f\l - x)- a ^ e+1 ^dx < C, 
1 p- 39 - 1 ( f'^dyuldy) dx 



and 



o 



< C 



J v 

1 p 1+9 (d x u) 2 dx^j 2 J\l - x)~ Q(1+3e) Qf \l - y) a{e ~ l) dy^ 1 dx 

< C Qf * p 1+9 (d x u) 2 dx^j J\l- x) 2 ^ 3+9 Ux 

< C (J p 1+9 {d x u) 2 dx^j , when 9 < 1, 



l / r i n 4 

35-1 / / I „0-l„ 



P [J \P u\dy) dx 

< c(j\ 9 -Wdx\ 2 j\l-x)-< 1+ ^ 

< C U 1 p^-Vdx) 2 j\l - x) 2 - a( ~ 3+9 Ux 



< C[ I P °- 1 u 2 dx ) , when < 1, 



1 p- 39 - 1 



(Jjl-yTV-Vdy) dx 



J \p 9 d y u\ + \p 9 1 u\dy^j dx 
< C ^ p 1+9 {d x u) 2 dx^j + C ( K J o A Vcfa^ , when 9 > 1. 
Thus, we get (|5.9p . Similarly, using a < ^ and a < 2 (i+g) » we can obtain (|5,10p easily. □ 
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Now, we can prove Lemma 12.51 as follows: 
Proof. Differentiating (|1.9j) -? respect to t, we get 
dttu 

= d t (r"" 1 ((2d + c 2 )/+V~Ms " P 1 ) J " 2d(n - 1)$ (r n - 2 ud x p 6 ) + 

:= I + J, (5.12) 

where I = r n ~ 1 d x ((2d + 02)p fl+1 § B (r n ~ 1 $u)) - 2ci(n - 1) (r n - 2 d t ud x p e ) . 
Multiplying (|5.12p by dtu and integrating it over x from to 1, we have 

4 f l(d t u) 2 dx + (-c 1 + c 2 ) I P e+l {d x {r n - l d t u)) 2 dx 
at ./o 2 n J n 



+2 (n l l Cl f p e+l {r n - l d tx u- p^r^dtuf dx = [ Jd t udx. 

n Jo Jo 



Denote 

B\ = p e +\^-\d tx u)\ B 2 = p e -\-\d t uf. 

Thus 

' ' ' o e+3 {d x u) 4 r 4n - 4 
-p e+ \d x u)W n - 4 + p^\d x u)\^ + p 2 (d x u) 2 



J Jd t udx < r] j (B 2 + B 2 )dx + C v J (p 6+3 {d x u) 4 r 4 
+p e+ \d x u) 2 u 2 r 2n - 4 + p 2 ^ e+l {d x u) 2 r 2n - 
+p 3 - e (d x u) 2 r 2n - 2 +p e - 1 u 4 ^j dx + J \d t ud t f\dx. 

Using a < r < R, m > , 2a(l — 9) > —1, Lemmas 12.21 and l5T2l we get 

J 1 p e+ \d x u) 4 r 4n - 4 dx < C (£ ( P 1+e {d x u) 2 + p°- l u 2 + {d t u) 2 ) dx\ 2 + C, 

I p e+ \d x u) 2 u 2 r 2n - 4 dx<C f p e+ \d x u) 4 dx + C f p e - l u 4 dx 
Jo Jo Jo 

< C f p e+3 {d x u) 4 dx + C [ pV-V^dx + C [ u 4m dx 

Jo Jo Jo 

< C^( K p l+e {d x u) 2 + p e - l u 2 + (d t u) 2 )dx^j + C, 

f 1 p 2 ~<- e+ \d x u) 2 r 2n ~ 2 dx < C f 1 p e+3 {d x u) 4 dx + C f 1 p^-^dx 
Jo Jo Jo 

< C^( y p l+e {d x u) 2 + p e - 1 u 2 + {d t u) 2 )dx^j +C, 



l f i r i 

p 3 ~ e (d x u) 2 r 2n - 2 dx <C p 4 (d x u) 4 dx + C p 2 - 29 dx 
o Jo Jo 
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< C (T (p 1+9 {d x u? + /-V + {d t uf) dx^j + c, 



< 



and 



p 2 (d x u) 2 dx <C f p e+3 (d x u) 4 dx + C [ p x - 6 dx 
o Jo Jo 

C U 1 (p 1+e {d x uf + A V + (dtu) 2 ) dx^j 2 + C, 



p e - 1 u i dx <C p^ 1] ^dx + C u Am dx < C. 



Since £ {r 11 ' 1 d x [(2 Cl + c 2 )p e+1 d x (r n - 1 u) - p~<} - 2c x (n - l)r n - 2 ud x p e + f} 2 dx\ t=0 is bounded 
(see (A3)), we have Jq u 2 (x, 0)dx is bounded. Thus, we have 



t pi r 



(dtu) dx + 

JO JO 

t / r l 

<c + c 



6+1 2n-2/o \2 , 8-1 -2/o \2 

2 



dxds 



f (J (p 1+9 (d x u) 2 + A V + (9^) 2 ) dx) A. 



From (|5.4p and (|5.9p . we have 

/ 

•70 



t r l 



+ 



JO 



■t / rl 



dxds 



< 



Cu + CkJ (J {p l+9 ^ u ? + P 9 ' 1 ™ 2 + (dtu) 2 )dx^j ds. 



Using Gronwall's inequality, we get 

■ i 



o 



(p^ 6 {d x uf + p 6 -\ 2 + {d t uf)dx 



t rl 



+ 



p" +l r 2n - 2 {d tx uY+p"-'r-- 2 (d t u) 



l„-2/ 



< 



10 Jo 

1 



dxds 



C24 



where 



< 2C 24 := C 6 , < t < T 4 , 



1 



2C24C25 

Using Lemmas I5.1N5.2I and (|5.13p . we can get (|2.15|W(|2.16p immediately. 



(5.13) 

(5.14) 
□ 



5.5 Proof of Lemma 12.61 

Proof. From (|5.1ip . we have 



\d x u\ Xo <c(p x ^- e -^ + \p- l u\ Xo + 



-0-1 



\dtu\dx 
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Ao 




Ao 




Ao\ 




+ 


p- - 1 f \p°d x u\dx 


+ 


p- 6 - 1 f\up e - l \dx 




J X 




J X 




J X 





for some positive constant C. 

From (fTTSl) . (l2i)l) . ([2~T9i Lemmas O and EE we have 

/ Ipl^^-^ds < C /* (1 - a;)^^- 9 - 1 ^ < C, 
Jo Jo 



\p- l u\ Xo dx< (J p 4 ™->o Xo dx 



1 X ^ 

1 \ 4m 

4m, 



u* m dx < a 



\d t u\dy 



Ao 



dx< f [p- - 1 ^- x^l^dxWdtuW^l 



< C / (1 - sjC-^+^+i^dx < C, 

■/ 







Ao ,-l 






p- 9 - 1 ! \p e d x u\dy 


dx< 


p- e -\ [ 


/o 


J X 


Jo 


Jo 



p- 8 -\ P e+ \d x ufdy)\{ p^dy) 



Ao 



dx 



<C (1 - x)- Xo ^ e+1 Ux < C, 
Jo 



and 







Ao 










p- 9 - 1 f \p d - l u\dy 


dx < 


f 


P- 1 


'-Hf'p 


10 


J X 




Jo 




Jo 



Ao 



< C / (1 - x)- Xoa{e+1) dx < C, 



where we use the fact that Ao < mrn { 40^+1 ; a (o+i) }• 

From above all, we can find a positive constant Cs such that Jq \d x u\ x °dx < C$. Using 
Lemma [2T2| (|2. 18j) and Sobolev's embedding theorem, we can get (|2.17p immediately. □ 
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